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Abstract. It is argued that accretion discs in young stellar objects may have hot coronae 
that are heated by magnetic reconnection. This is a consequence of the magneto-rotational 
instability driving turbulence in the disc. Magnetic reconnection away from the midplane 
leads to heating of the corona which, in turn, contributes to driving disc winds. 
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1. Introduction 

Accretion discs provide the channel through 
which matter from the outside is able to col- 
lapse onto some common central object such 
as a protostar. Prior to the collapse, most of the 
matter will be in random motion and will have 
excess angular momentum with respect to the 
forming star The disc allows the matter to get 
rid of this excess angular momentum, but in or- 
der to do this there needs to be some enhanced 
viscosity of some type. It is now generally 
accepted that this is provided by turbulence 
which is driven by the magneto-rotational in- 
stability (MRI; see Balbus & Hawley 1991, 
1998). 

An automatic by-product of MRI-driven 
turbulence is the possibility of significant re- 
sistive Joule heating in the layers aways from 
the midplane, where the density is small and 
hence the heating per unit mass is large. This is 
demonstrated in Fig.[Tl where we show that the 
total (magnetic and kinetic) stress, that leads 
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to viscous and resistive heating, does not de- 
crease with height away from the midplane, 
even though the product of rate of strain and 
density does decrease because of decreasing 
density. This was originally demonstrated only 
for nearly isothermal discs (Brandenburg et al. 
1996a), see Fig. [T] but this has now also been 
shown for discs with full radiation transfer in- 
cluded (Turner 2004). 

An important expectation from this is that 
the discs should be surrounded by hot coronae. 
In other words, not only the central protostar, 
but also the disc itself should be surrounded by 
hot coronae that are heated by the associated 
Joule (or resistive) heating via magnetic recon- 
nection, as has been demonstrated in numeri- 
cal experiments (Galsgaard & Nordlund 1996). 
This leads to the possibility of hot coronae with 
temperatures close to the Virial temperature, 
Tvir, given by 

CpT^Ar) = GMJr, (1) 

where Cp is the specific heat at constant pres- 
sure, G is Newton's constant, M, is the mass 
of the protostar, and r is the distance from the 
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Fig. 1. Dependence of tlie stress component f^^ 
(here denoted by Ttot), separately for the kinetic and 
magnetic contributions, t^^^ and T„,ag, respectively, 
together with the sum of the two denoted by total 
(left) as well as the vertical dependence of density 
and sound speed (right). Note that Ttot is neither pro- 
portional to the density p nor the sound speed c^. 
[Adapted from Brandenburg et al. (1996a).] 



central object. The only process that can pos- 
sibly counteract the tendency for strong coro- 
nal heating is radiative cooling. At present it is 
unclear what temperatures will be achieved if 
realistic cooling is included. 

2. The mean accretion stress 

The basic dynamics of the MRI can be mod- 
eled using simulations in a local cartesian ge- 
ometry using shearing sheet boundary condi- 
tions in the radial direction. One of the main 
output parameters is the Shakura-Sunyaev 
viscosity parameter, ass, which is a non- 
dimensional measure of the turbulent viscosity 
Vt, in terms of the sound speed Cs and the disc 
scale height H, i.e. 



Vt = assCsH. 



(2) 



Assuming that there is no external torque, this 
ass is given by 



-pVtW 



do. 

dm' 



(3) 



where m is the cylindrical radius, p is the den- 
sity, Q is the mean angular velocity, and f^if, is 
the total 'horizontal' stress given by 



(4) 



where b and u denote the fluctuating compo- 
nents of the magnetic and velocity fields, and 
fiQ is the vacuum permeability. The first term. 
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Fig. 2. Time series of the Shakura-Sunyaev viscos- 
ity alpha, cFss- Here, time is normalized in terms of 
the orbital time, Ty^t = 2n/Sl, where Q is the local 
angular velocity. Note that ctss fluctuates strongly 
in time about an average value of around 0.01. 
[Adapted from Brandenburg (1998).] 



i.e. the magnetic stress, is usually the largest 
term. The resulting stress determined from the 
simulations and expressed in nondimensional 
form as ass(t), is shown in Fig. [2] The simula- 
tion domain covers the upper disc plane only. 

The expression in equation (O would be 
used in the evolution equation for the mean an- 
gular velocity, Q, where the term Vt would sim- 
ply add to the microscopic viscosity v to give a 
total viscosity, vj = Vt + v. On the other hand, 
equation (HJl would be used to determine t^^, 
and hence ass, from three-dimensional simu- 
lation data. 

The value of ass(0 is notoriously small 
(around 0.01; see Hawley et al. 1996, Stone 
et al. 1996, Brandenburg et al. 1996b). 
Observational evidence suggests a typical 
range ass ~ 0.1. ..0.4 (King et al. 2007). Larger 
value of ass can be obtained if there is an ex- 
ternally imposed magnetic field. Also global 
simulations tend to yield larger values (Hawley 
2000). 

The variation of ass(0 shows a typical time 
scale of around 15 orbits. This is here a conse- 
quence of the emergence of cyclic large scale 
dynamo action. The magnetic field reverses be- 
tween two maxima, so the actual period is ac- 
tually around 30 orbits. Physically, this time 
is related to the turbulent diffusion time. This 
is now reasonably well understood and the 
details of the dynamo depend sensitively on 
the boundary conditions. Such cycles are not 
expected in global geometry (Brandenburg 8c 
Campbell 1997). 
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3. The mean electromotive force 
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An important goal is to simulate outflows from 
discs that can produce a magnetic field. In 
other words, the magnetic field necessary for 
launching outflows is now no longer assumed 
as being given, but it is self-consistently gener- 
ated in the disc, albeit in parameterized form. 
This feature was investigated by von Rekowski 
et al. (2003), who studied the possibility of col- 
limated outflows in the absence of an ambi- 
ent magnetic field. This idea seems now quite 
reasonable also from an observational point of 
view. Menard & Duchene (2004) report obser- 
vations of the Taurus-Auriga molecular cloud 
which is a nearby star-forming region of low- 
mass stars. They find that jets and outflows 
from classical T Tauri star-disc systems in 
this cloud are not always aligned with the lo- 
cal magnetic field of the cloud where they are 
born. This suggests that jets are not likely to be 
driven by the cloud's own field. However, theo- 
retical models have shown that magnetic fields 
are important for jet acceleration and collima- 
tion. Furthermore, there is now also observa- 
tional evidence that jets and outflows are essen- 
tially hydromagnetic in nature. Observations 
of CO outflows exclude purely radiative and 
thermal driving of jets and outflows in regions 
of star formation of low-mass stars (Pudritz 
2004). Thermal pressure might still be dynam- 
ically important to lift up the winds to the criti- 
cal point. Moreover, there are observational in- 
dications of warm wind regions with tempera- 
tures of a few times 10"^ K (Gomez de Castro 
& Verdugo 2003). These winds cannot be due 
to stellar radiation because the star is too cool, 
suggesting that the warm winds are driven 
magnetohydrodynamically. So it is quite plau- 
sible that a disc dynamo might be necessary for 
the launching of jets and outflows. 

In order to address problems such as the 
launching of disc winds from hot disc coronae, 
it is preferable to solve the equations in global 
geometry. This increases the computational de- 
mands, but some progress can be made by re- 
stricting oneself to axisymmetry. On the other 
hand, dynamo action is impossible in axisym- 
metry (Cowling 1933). However, it is sufficient 
to capture the essential three-dimensional ef- 
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Fig. 3. Scatter plot of the volume averaged toroidal 
electromotive force, (fi^), versus the volume aver- 
aged toroidal magnetic field strength {B^). Here, 
(G^) is normalized with respect to the rms veloc- 
ity "V and the equipartition field strengths B^q = 
^j^Qpu-. Both {&^) and {B^) are obtained at different 
times during the simulation. Here yUo is the vacuum 
permeability. Note the negative slope in the diagram 
around {B^) = 0, indicating that the a effect in mean 
field dynamo theory is negative. The solid line rep- 
resents a fit to a standard quenching law of the form 
-/3/(l + fS^), where yS = (B^)IB^^. [Adapted from 
Brandenburg & Donner (1997).] 



fects in a statistical sense by considering the 
averaged equations. 

In the averaged equations the nonlinear 
terms lead to correlations, just like the stress 
terms considered in equation (|4]i. In the present 
case there is the mean electromotive force. 



S-uxb, 



(5) 



where lower case u and b denote fluctuations, 
i.e. velocity and magnetic fields are split into 
mean and fluctuating components via U = U + 
u and B - B + b. Equation (|5]l is used to deter- 
mine £ from simulation data. This expression 
is thus analogous to equation (|4]i. 

In isotropic turbulence £ = (just like 
there is no off-diagonal terms in the stress con- 
sidered above). However, in the presence of ro- 
tation and stratification £ can be finite (Krause 
& Radler 1980). In particular, theory suggests 
that £ has a component along the direction of 
the mean field, i.e. 



£ - QdynB + 



(6) 
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where adyn is the dynamo alpha which has the 
dimensions of a speed. The dots indicate the 
presence of additional terms such as turbulent 
resistivity, for example. The expression equa- 
tion (|6]l would be used in a mean field model 
in the averaged induction equation for B, and 
is thus analogous to equation (O which is used 
in the averaged momentum equation. 

In Fig. [3] we show a scatter plot of 
the toroidal electromotive force, (£0), versus 
the volume averaged toroidal magnetic field 
strength {B^} that has been obtained from a 
simulation of MRI-driven turbulence in the up- 
per disc plane. Note the negative slope in the 
diagram indicating that the a effect in mean 
field dynamo theory is negative (Brandenburg 
et al. 1995, Ziegler & Rudiger 2000). This is an 
unexpected result, because theory predicts ffdyn 
to be positive in the upper disc plane. Various 
proposals for this negative sign have been put 
forward (Brandenburg 1998, Rudiger & Pipin 
2000, Blackman & Tan 2004). More detailed 
analysis confirms the negative value of ffdyn 
and is also able to give its dependence on the 
distance from the midplane and other aspects 
(Brandenburg & Sokoloff 2002). 

There is a great deal of work concerned 
with the possibility of so-called catastrophic 
quenching of adyn that depends on the value 
of the microscopic magnetic Reynolds number. 
We cannot go into these aspects here, but we 
mention only that such quenching can be al- 
leviated by magnetic helicity fluxes out of the 
domain (Blackman & Field 2000, Kleeorin et 
al. 2000), and that such fluxes can be driven 
by shear (Vishniac & Cho 2001, Subramanian 
& Brandenburg 2004, 2006). The end re- 
sult might be an a effect that is no longer 
catastrophically quenched; see Brandenburg & 
Subramanian (2005) for a review. In the fol- 
lowing we will just use this as an established 
fact, as has been done elsewhere in the litera- 
ture (Campbell 2000, Bardou et al. 2001). 

4. Structured outflows from discs 

A mean field dynamo model with piece-wise 
polytropic hydrodynamics was proposed by 
von Rekowski et al. (2003). Like Ouyed & 
Pudritz (1997a,b) they start with an equilib- 
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Fig. 4. General structure of the outflows typi- 
cally obtained in the model of von Rekowski 
et al. (2003), where the cool, dense disc 
emits a thermally driven wind and a magneto- 
centrifugally driven outflow in a conical shell. 
[Adapted from von Rekowski et al. (2003).] 



rium corona, where they assume constant en- 
tropy and hydrostatic equilibrium according to 
CpT(r) - GMt/r. In order to make the disc 
cooler, a geometrical region for the disc is pre- 
scribed (see Fig. |4|. An entropy contrast be- 
tween disc and corona is chosen such that the 
initial disc temperature is about 3 x lO-^K in 
the bulk of the disc. The low disc tempera- 
ture corresponds to a high disc density of about 
10"'° . . . 10"^ g cm"^. For the disc dynamo, the 
most important parameter is the dynamo coef- 
ficient ffdyn in the mean field induction equa- 
tion. The a effect is antisymmetric about the 
midplane and restricted to the disc. 

The upper panel of Fig.|5]illustrates that an 
outflow develops that has a well-pronounced 
structure. Within a conical shell originating 
from the inner edge of the disc, the terminal 
outflow speed exceeds 500 km/s, and tempera- 
ture and density are lower than elsewhere. The 
inner cone around the axis is the hottest and 
densest region where the stellar wind speed 
reaches about 150 km/s. The wind that devel- 
ops from the outer parts of the disc has inter- 
mediate values of the speed. 

The structured outflow is driven by a com- 
bination of different processes. A significant 
amount of angular momentum is transported 
outwards from the disc into the wind along 
the magnetic field, especially along the strong 
lines within the conical shell (see the lower 
panel of Fig. |5]l. The magnetic field geom- 
etry is such that the angle between the ro- 
tation axis and the field lines threading the 
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Fig. 5. Upper panel: poloidal velocity vec- 
tors and poloidal magnetic field lines (white) 
superimposed on a color scale representa- 
tion of log/i. Specific enthalpy h is di- 
rectly proportional to temperature T, and 
log/i - (-2,-1,0,1) corresponds to T ^ 
(3xl0^3xlO'*,3xl0^3xlO^)K. The black 
dashed line shows the fast magnetosonic sur- 
face. The disc boundary is shown with a thin 
black line, the stellar surface is marked in 
red. The dynamo a^yn coefficient is negative in 
the upper disc half, resulting in roughly dipo- 
lar magnetic symmetry. Averaged over times 
t ^ 897 . . . 906 days. Lower panel: color scale 
representation of the specific angular momen- 
tum, normalized by the maximum angular mo- 
mentum in the disc, with poloidal magnetic 
field lines superimposed (white). The black 
solid line shows the Alfven surface, the white 
dashed line the sonic surface. Same model as in 
the upper panel and averaged over same times. 
[Adapted from von Rekowski & Brandenburg 
(2004).] 



disc exceeds 30° at the disc surface, which 
is favorable for magneto-centrifugal accelera- 
tion (Blandford & Payne 1982). However, the 
Alfven surface is so close to the disc surface 
at the outer parts of the disc that acceleration 
there is mainly due to the gas-pressure gradi- 
ent. In the conical shell, however, the outflow 
is highly supersonic and yet sub-Alfvenic, with 
the Alfven radius a few times larger than the 
radius at the footpoint of the field lines at the 
disc surface. The lever arm of about 3 is suf- 
ficient for magneto-centrifugal acceleration to 
dominate in the conical shell (cf. Krasnopolsky 
et al. 1999). 

5. Star-disc coupling 

The interaction of a stellar magnetic field with 
a circumstellar accretion disc and its magnetic 
field was originally studied in connection with 
accretion discs around neutron stars (Ghosh 
& Lamb 1979), but it was later also applied 
to protostellar magnetospheres (Konigl 1991, 
Cameron & Campbell 1993, Shu et al. 1994). 
Most of the work is based on the assumption 
that the field in the disc is constantly being 
dragged into the inner parts of the disc from 
large radii. The idea behind this is that a mag- 
netized molecular cloud collapses, in which 
case the field in the central star and that in 
the disc are aligned (Shu et al. 1994). This 
was studied numerically by Hirose et al. ( 1 997) 
and Miller & Stone (1997). In the configura- 
tions they considered, there is an X-point in the 
equatorial plane (see left hand panel of Fig. |6]l, 
which can lead to a strong funnel flow. 

The other alternative has been explored by 
Lovelace et al. (1995) where the magnetic field 
of the star has been flipped and is now anti- 
parallel with the field in the disc, so that the 
field in the equatorial plane points in the same 
direction and has no X-point. However, current 
sheets develop above and below the disc plane 
(see right hand panel of Fig. |6]l. This configu- 
ration is also referred to as the X-wind model. 
Ironically, this is the field configuration with- 
out an X-point. 

Numerical simulations of such a field con- 
figuration by Hayashi et al. (1996) confirm the 
idea by Lovelace et al. (1995) that closed mag- 
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Fig. 6. Sketcii siiowing the formation of an X-point 
when the disc field is aligned with the dipole (on 
the left) and the formation of current sheets with no 
X-point if they are anti-aligned (on the right). The 
two current sheets are shown as thick lines. In the 
present paper, the second of the two configurations 
emerges in all our models, i.e. with current sheets 
and no X-point. [Adapted from von Rekowski & 
Brandenburg (2004).] 



netic loops connecting the star and the disc 
are twisted by differential rotation between the 
star and the disc, and then inflate to form open 
stellar and disc field lines (see also Bardou 
1999, Agapitou & Papaloizou 2000, Uzdensky 
et al. 2002). Goodson et al. (1997,1999) and 
Goodson & Winglee (1999) find that for suffi- 
ciently low resistivity, an accretion process de- 
velops that is unsteady and proceeds in an os- 
cillatory fashion. The same result has also been 
obtained by von Rekowski & Brandenburg 
(2004). It should be emphasized that in their 
case the resulting field geometry is always the 
second one in Fig. |6] i.e. the one with cur- 
rent sheets and no X-point. The inflating mag- 
netosphere expands to larger radii where mat- 
ter can be loaded onto the field lines and be 
ejected as stellar and disc winds. Reconnection 
of magnetic field lines allows matter to flow 
along them and accrete onto the protostar, in 
the form of a funnel flow (see also Romanova 
et al. 2002). Consequently, the outflows show 
episodic behavior; see also Matt et al. (2002) 
and von Rekowski & Brandenburg (2004). 

In the model of von Rekowski & 
Brandenburg (2004) the highly episodic accre- 
tion flow and accretion rate are correlated with 
the configuration of the magnetosphere. Low 



values of at the position (m,z) - (0.6,0.4), 
correspond to a configuration were closed 
magnetospheric lines penetrate the inner disc 
edge, thus connecting the star to the disc (as in 
the left hand panel of Fig.|7]i. In this case, disc 
matter is loaded onto magnetospheric field 
lines and flows along them to accrete onto the 
star As a consequence, the accretion rate is 
highest during these phases, typically up to 
between lO^^^Moyr"' and 2.5 x 10"**Moyr"'. 

High values of B~ at the position (m, z) - 
(0.6, 0.4), coiTespond to a configuration when 
the outer field lines of the magnetosphere have 
opened up into disconnected open stellar and 
disc field lines, thus disconnecting the star 
from the disc (as in the right hand panel of 
Fig- El- In this case, matter is lost directly into 
the outflow and there is no net accretion of disc 
matter. 

6. Conclusions 

Simulations of disc dynamos as well as of 
star-disc coupling have shown several surprises 
that were not originally expected. On the other 
hand, some of these results are still only ten- 
tative and one cannot really be sure to what 
extent they are artifacts of the model setup or 
other technical limitations. The fact that adyn 
is negative in the upper disc plane is one such 
result that has so far only been obtained in 
the local shearing box simulations. It would 
be good to check this result independently us- 
ing global models. Most of the recent global 
simulations tend to focus on the magnetic field 
structures and their effects on the gas rather 
than on studying the magnetic field generation, 
for example by identifying parameterizations 
that could be used for mean field modeling. 

Regarding star-disc coupling there are 
questions concerning the strength of numeri- 
cal viscosity and other technical aspects that 
make it difficult to be sure about the signifi- 
cance of the result. Here, however, the basic 
processes involved are now well understood, 
in particular the inflation of field lines that 
are being differentially sheared are they per- 
meate the disc. As argued by Matt & Pudritz 
(2004, 2005), this makes the disc-locking of 
the star rather inefficient. This motivates the 
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Fig. 7. Colour/grey scale representation of tlie density (bright colors or light shades indicate high values; 
dark colors or dark shades indicate low values) with poloidal magnetic field lines superimposed (white) and 
the azimuthally integrated mass flux density, represented as the vector In-mgiu^, ii-), shown with arrows 
(except in the disc where the density is high and the mass flux vectors would be too long). The left hand 
panel is at the time t = 9S when the star is magnetically connected to the disc and the right hand panel is at 
the time t = 102 when the star is magnetically disconnected from the disc. [Adapted from von Rekowski & 
Brandenburg (2004).] 



search for other mechanisms facilitating stel- 
lar braking. A plausible candidate is braking 
by a stellar wind, which was also found by von 
Rekowski & Brandenburg (2006). Extending 
these types of calculations to three dimensions 
and to higher resolution remain important fu- 
ture tasks (see, e.g., von Rekowski & Piskunov 
2006 for initial studies). 
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